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The  orbit  of  China  6  rocket,  1976-87B,  has  been  determined  at  51  epochs 
during  its  17-month  life,  using  the  RAE  orbit  refinement  computer  program,  PROP  6, 
with  over  4000  radar  and  optical  observations  from  49  stations.  The  orbital 
accuracy  is  about  100  metres,  radial  and  cross- track,  on  average. 

The  rotation  rate  of  the  upper  atmosphere,  A  rev/day  ,  for  the  height-band 
200-230  km,  was  calculated  from  the  decrease  in  orbital  inclination  (after  being 
cleared  of  perturbations)  to  give  the  following  results: 

(1)  for  morning  conditions,  A  =  0,9  for  May- June  and  August-September  1977, 
at  215  km  mean  height;  A  =  0.7  for  October-November  1977,  at  210  km; 

A  =  0.8  ±  0.05  for  January-February  1978,  at  200  km; 

(2)  for  evening  conditions,  A  =  1.2  for  July  and  September-October  1977,  at 
215  km; 

(3)  for  mean  (morning  plus  evening)  conditions,  A  =  1.0  ±  0.1  between  October 
1976  and  May  1977,  at  230  km;  A  -  0.8  ±  0.1  for  December  1977-January 
1978,  at  215  km  and  mean  latitude  57°S. 

Values  of  density  scale  height  have  been  obtained  from  the  variation  in 
perigee  height,  including  several  values  during  the  final  16  days  to  decay. 
Comparison  with  CIRA  1972  values  shows  agreement  mostly  within  10  per  cent. 
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1  INTRODUCTION 

China  6  rocket  (1976-87B)  was  launched  on  30  August  1976,  and  remained  in 
orbit  for  about  17  months  before  decaying  in  the  Earth's  atmosphere  on 
4  February  1978.  The  initial  orbit  was:  inclination  69.2  degrees,  perigee 
height  190  km,  apogee  height  2150  km,  eccentricity  0.129  and  period  108.7  minutes. 
The  orbit  has  been  determined  at  35  epochs  from  more  than  1700  radar  and  optical 
observations  with  the  RAE  orbit  refinement  program* ,  in  the  PROP  6  version.  The 
highly— accurate  Hewitt  camera  observations  from  Malvern  were  used  in  six  of 
the  orbital  determinations.  Further  orbits  were  determined  for  each  o^  the  last 
16  days  of  the  life,  using  about  1300  observations  obtained  from  the  North 
American  Air  Defense  Command  (NORAD) ,  about  50  US  Navy  observations  and  a 
Hewitt  camera  run  of  observations. 

As  the  orbit  contracted  under  the  influence  of  air  drag,  it  passed  through 
14th  and  15th-order  resonances,  but  too  quickly  for  any  harmonic  coefficients 
to  be  determined  with  worthwhile  accuracy.  The  main  purpose  here  is  to  find 
the  upper-atmosphere  rotation  rate  from  changes  in  inclination;  the  resonant 
perturbations  are  required  to  improve  the  accuracy  of  the  rotation  rate,  which 
has  been  determined  before,  between  and  after  the  resonances. 

Upper-atmosphere  scale  heights  have  been  calculated  for  the  height-band 
200-230  km  and  compared  with  CIRA  1972  values. 

2  ORBIT  DETERMINATION 
2. 1  Observations 

Over  3600  radar  and  optical  observations  were  used  to  determine  the  orbit 
at  51  epochs  between  16  October  1976  and  4  February  1978.  A  further  400 
observations  of  low  elevation  were  not  used.  During  the  final  orbit  determina¬ 
tions,  about  700  observations  were  rejected,  leaving  about  3000  observations 
or  near  58  per  orbit,  on  average.  These  included  over  1100  from  the  US  Navy 
Navspasur  system  and  about  1400  from  the  assigned  and  contributing  sensors  of 
the  North  American  Air  Defense  Command  (N0RAD)  Space  Detection  and  Tracking 
System  (SPADATS).  A  further  300  were  visual,  from  volunteer  observers  reporting 
to  the  Appleton  Laboratory  at  Slough;  200  or  so  were  from  British  radar;  18  from 
the  Malvern  Hewitt  camera,  distributed  over  seven  orbits;  24  from  Jokioinen, 
Finland  and  15  from  the  kine theodolite  at  the  South  African  Astronomical 
Observatory. 


2.2  Observational  accurac- 


Table  1 

Residuals  for  selected  stations 


Station 

Number 

of 

observations 

Rms  residuals 

Range 

km 

Minutes  of 

arc 

RA 

r 

Dec 

Total 

1 

US  Navy 

84 

1.8 

2.4 

2.9 

2 

US  Navy 

33 

2.3 

2.1 

3.2 

3 

US  Navy 

42 

2.0 

1.7 

2.6 

4 

US  Navy 

51 

1.6 

2.1 

2.7 

5 

US  Navy 

95 

2.0 

2.2 

3.0 

6 

US  Navy 

1 10 

1.8 

2.4 

3.0 

29 

US  Navy 

693 

0.8* 

0.5* 

0.5* 

414 

Capetown 

10 

2.3 

3.0 

3.8 

2122 

Malvern  5 

27 

1.7 

1.9 

2.5 

2155 

Bahrein  2 

17 

1.6 

3.1 

3.5 

2265 

Farnham 

29 

2.0 

2.2 

KITH 

2303 

Malve-n  Hewitt  camera 

18 

0.08 

0.09 

2414 

Bournemouth 

103 

3.2 

3.8 

4.9 

2420 

Willowbrae 

1 19 

2.1 

2.1 

3.0 

2577 

Cape  kinetheodolite 

15 

0.6 

0.5 

0.8 

6702 

Jokioinen 

24 

3.5 

2.7 

4.4 

*  Geocentric 


Table  1  lists  the  observing  stations  which  have  contributed  at  least  10 

observations  to  the  determination  of  the  first  35  orbits.  The  observational 

accuracy  of  these  stations  is  given  in  the  form  of  rms  residuals,  obtained  from 
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the  RAE  computer  program  ORES  .  The  Hewitt  camera  observations  have  the  highest 
accuracy,  with  about  5  seconds  of  arc;  the  US  Navy  stations  have  topocentric 
accuracies  near  3.0  minutes  of  arc;  the  Jokioinen  theodolite  4.4  minutes  of  arc; 
the  Cape  kinetheodolite,  0.8  minutes  of  arc.  The  stations  with  visual  observers 
reporting  to  the  Appleton  Laboratory  at  Slough,  have  residuals  between  2.5  and 
4.9  minutes  of  arc.  Each  observer  has  been  sent  a  list  of  his  own  residuals. 

2.3  Orbital  accuracy 

The  computed  sets  of  orbital  elements  at  51  epochs  between  October  1976  and 
February  1978  are  shown  in  Table  2,  with  standard  deviations  where  appropriate. 
The  sd  in  inclination  varies  from  0.0005  to  0.0028°,  with  an  rms  of  0.0015°. 

For  the  final  16  daily  orbits,  the  sd  varies  from  0.0005  to  0.0010°,  with  an  rms 
of  0.0007°.  For  eccentricity,  plotted  in  Fig  1,  the  rms  for  the  51  values  of  sd 
is  2  x  10  ,  but  only  1  x  10  ^  (equivalent  to  70  m  in  perigee  height)  for  the 
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final  16  values.  For  argument  of  perigee,  w  ,  or  mean  anomaly  at  epoch,  , 
the  (similar)  sd  varies  from  0.005  to  0.090°.  The  values  of  orbital  decay 
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rate,  M2  ,  plotted  in  Fig  2,  have  sd  varying  from  0.0002  to  0.003  deg/day 
for  the  first  35  orbits;  but  for  the  final  16  orbits  (approaching  decay), 
shown  in  the  inset  on  Fig  2,  the  sd  increases  and  on  the  actual  decay 
day,  M2  *  30.71  ±0.05  deg/day^. 

3  VARIATION  IN  PERIGEE  HEIGHT 

It  is  useful  to  analyse  the  perigee  height  first,  because  the  results  are 
used  in  deciding  the  height  at  which  the  rotation  rate  applies  (see  section  6) . 

The  values  of  semi  major  axis,  a  ,  and  eccentricity,  e  ,  given  in 
Table  2,  were  used  to  calculate  perigee  height  over  a  spherical  Earth,  h^  , 
from 

h  =  a ( 1  -  e)  -  R  ,  (1) 

P 

where  R  is  the  Earth's  equatorial  radius,  taken  as  6378.14  km.  The  values, 
cleared  of  perturbations,  are  plotted  in  Fig  3,  up  to  MJD  43524,  as  crosses. 

The  main  oscillation  is  due  to  the  odd-harmonic  perturbation  in  e  ,  which, 

3 

combined  with  the  lunisolar  perturbation,  was  calculated  using  the  PROD  program  , 
to  give  a  total  perturbation  Ae  .  A  parameter  Q  ,  the  perigee  height  over  a 
spherical  Earth,  cleared  of  perturbations,  can  now  be  determined  from 


Q 


h  +  aAe  , 
P 


and  the  values  obtained  are  shown  in  Fig  3  by  triangles.  The  perturbation  Ae 
is  assumed  to  be  zero  initially,  so  that  the  Q  values  are  relative  (rather  than 
absolute).  This  is  not  significant  here,  since  only  the  slope  dQ/dt  is  used 
(in  section  4)  for  calculating  density  scale  height. 


The  values  of  h^  and  Q  in  the  final  16  days  before  satellite  decay  are 
plotted  in  Fig  4  (which  is  an  extension  of  Fig  3) . 


DENSITY  SCALE  HEIGHT 


The  density  scale  height  H  is  a  measure  of  the  rate  of  decrease  of  density 

p  as  height  y  increases,  and  is  defined  by  i  •  Values  of  H  were 

.  H  P  dy  4 

calculated  from  Q  ,  the  rate  of  change  of  Q  due  to  air  drag,  where  ,  for 
ae/H  >  3  , 
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.  2  \ 
sin  i  cos  2o)] 


e'  is  the  ellipticity  of  the  atmosphere  (*0.00335)  and  H  is  at  a  height 
1.5  above  the  satellite's  perigee  height  y^  .  The  other  variables  are 
defined  in  Table  2.  Ignoring  small  perturbations,  we  may  take  y^  as  given  by 
the  right-hand  side  of  equation  (1)  with  R  replaced  by  the  local  Earth  radius 
at  perigee  latitude.  For  1976-87B,  at  69°  inclination. 


y  =  h  +  18.67  sin  w 
P  P 


Values  of  Q  for  use  in  equation  (2)  are  obtained  in  the  form  AQ/At  , 

from  the  change  AQ  in  Q  over  a  suitable  time  interval  At  .  Since  Q  has 

an  accuracy  of  about  0.2  km,  values  of  AQ  of  3  km  or  more  are  required  to  give 
« 

values  of  Q  accurate  to  5  per  cent,  and  the  time  intervals  were  chosen  on  this 
basis.  Values  of  H  were  obtained  from  equation  (2)  using  these  values  of  Q 
and  mean  values  of  the  other  parameters.  The  values  of  H  are  plotted  against 
time  span  in  Fig  5,  as  circles  with  sd,  up  to  MJD  43540.  At  this  point, 
ae/H  -  3. 

For  the  final  three  (daily)  calculations  of  H  ,  the  value  of  ae/H 
was  less  than  3.  The  'phase  2'  regime  discussed  in  Ref  5  (page  88)  can  then  be 
used:  equation  (5.35)  of  Ref  5  gives 

-  y0  +  le(4  -  3y^  -  yQy2)  -  ic  cos  2u)(yQ  -  2y2  +  yQy3)  =  3  ,  say  (3) 


e'a(l  -  e)  .  2 

whare  c  *  - ~ -  sin  i  ,  y  *  I  /I  ,  and  I  is  the  Bessel  function  of 

zn  r  r  J  r 

the  first  kind  and  imaginary  argument,  of  order  r  {ibid,  page  36). 


To  calculate  H  from  Q  ,  put 


q  =  ^2  -  da  = 

v  da  dt 


d(a  -  x)  ,  2na 
da  3n  ’ 


where  n(=Mj)  is  the  mean  motion  and  n  =  2M2  .  Hence 


rjj  =  l/(l  +  3MjQ/4M2a)  =  a  ,  say. 
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H  is  now  determined  by  selecting  two  likely  values,  and  ,  which  are 

used  to  calculate  gj  and  from  equation  (3),  and  then  by  linear  interpolation 
from 


H,(02  -  a)  +  H2(a  -  6,) 


These  three  values  of  H  are  included  in  Fig  5  as  circles,  with  sd,  for 
43540  ^  MJD  ^  43543.  The  sd  for  all  the  H  values  shown  have  been  estimated 
from  the  errors  in  perigee  height  (usually  0.2  km).  The  height  y^  at  which 
H  applies  has  been  plotted  in  Fig  5;  also  the  10.7  cm  solar  flux,  F^  y  ,  in 
the  form  of  mean  monthly  values  up  to  MJD  43529,  and  daily  thereafter. 

For  comparison,  values  of  scale  height  H  obtained^  from  CIRA  1972,  for 
the  same  heights  and  exospheric  temperatures  as  the  calculated  values,  have  also 
been  plotted  in  Fig  5  as  crosses.  These  observational  values  of  scale  height, 
between  180  and  240  km  height,  are  about  7  per  cent  rms  above  CIRA  1972  values  for 
high  solar  activity  in  1977-8.  In  a  previous  study^,  the  observational  values 
(at  similar  heights)  were  about  10  per  cent  below  CIRA  1972  values  for  lower 
solar  activity  (1972-6). 

The  observational  density  scale  height  H  has  been  replotted  against 
height  in  Fig  6,  with  comparative  curves  from  CIRA  1972,  for  various  exospheric 
temperatures.  The  calculated  exospheric  temperatures  are  shown  in  brackets 
after  the  observational  values  (indicated  by  crosses  generally  or  by  circles 
near  decay).  The  calculations  were  made  using  CIRA  1972  with  the  appropriate 
values  of  solar  10.7  cm  radiation  energy  and  geomagnetic  index,  with  adjustments 
for  semi-annual  variations.  Again  the  observational  values  of  H  are  higher 
than  the  values  given  by  CIRA  1972 ,  for  the  appropriate  exospheric  temperature , 
the  difference  being  6  per  cent  rms. 

5  RESONANCES 

5. 1  14th-order 

The  orbit  of  China  6  rocket  was  perturbed  by  I4th-order  resonance  between 
May  and  August  1977,  with  exact  14th-order  resonance  occurring  on  25  June 
(MJD  43319).  Unfortunately  there  were  only  five  PROP  values  of  inclination 
during  this  time-span  for  determining  lumped  14th-order  geopotential  coefficients; 
a  further  six  US  Navy  values  of  inclination  were  also  available  -  but  even  11 
values  are  insufficient  to  determine  such  coefficients  accurately.  However,  a 
fitting  was  made  to  these  11  values,  after  removing  unwanted  perturbations,  as 
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the  resonance  effect  on  inclination  is  required  for  improving  the  accuracy  in 

calculating  the  rotation  rate  of  the  atmosphere  (section  6).  The  perturbations 

3 

to  be  removed  were  found  using  the  computer  program  PROD  and  the  resulting 

g 

inclinations  were  fitted  using  the  computer  program  THROE  by  the  method 
described  in  Ref  7,  using  only  the  (y,q)  51  (1,0)  terms.  The  inclination  change 
due  to  14th-order  resonance  was  found  to  be  0.0058°,  and  is  shown  plotted  in 
Fig  7  as  a  vertical  line  at  exact  14th-order  resonance. 

5.2  15th-order 

The  change  in  inclination  at  29:2  resonance  is  too  small  to  be  worth 
considering,  but  at  15th-order  resonance  the  orbit  was  perturbed  sufficiently 
to  warrant  an  estimation  of  an  allowance  to  be  made  in  the  inclination  fitting 
of  Fig  7.  The  resonant  perturbation  is,  however,  too  rapid  for  the  determination 
of  lumped  geopotential  coefficients.  It  occurred  between  October  and  December 
1977,  with  exact  15th-order  resonance  on  25  November. 

There  were  only  seven  PROP  and  two  US  Navy  values  of  inclination  available 
and  the  fitting  was  made  using  (y,q)  =  (1,0)(1,±1).  The  overall  change  in 
inclination  due  to  15th-order  resonance  was  0.0022°  and  is  plotted  in  Fig  7  as 
a  vertical  line  at  exact  15th-order  resonance  (MJD  43472).  The  method  of 
fitting  is  as  described  in  Ref  7. 

6  ATMOSPHERIC  ROTATION  RATE 

The  51  derived  values  of  inclination,  given  in  Table  2,  were  cleared  of 

3 

lunisolar  and  geopotential  perturbations  using  the  PROD  program  with  numerical 

integration  at  one-day  intervals.  These  modified  values  are  plotted  in  Fig  7. 

The  change  in  inclination  was  then  calculated  for  several  values  of  atmospheric 

rotation  rate  (expressed  as  A  times  the  Earth's  rotation  rate,  in  rev/day) , 

9  ... 

using  oblate-atmosphere  theory  ,  with  numerical  integration  at  about  20-day 
intervals  (corresponding  to  22.5°  steps  in  argument  of  perigee,  to  ).  For  the 
final  16  days  of  the  orbit,  however,  the  intervals  were  reduced  to  1  day. 

For  the  first  7  months,  MJD  43067-43282,  the  local  time  at  perigee  (see 
scale  at  top  of  Fig  7)  passed  through  morning  and  evening  At  least  twice 
each,  so  the  value  of  A  obtained  is  an  average  one  and  is  given  by 
A  *  1.0  ±0.1  rev/day  for  a  mean  height  of  y(=  y^  +  0.75H)  =  228  km. 

In  the  last  16  days,  where  there  are  daily  orbits,  the  values  of  i  , 
which  have  been  replotted  with  greater  accuracy  in  Fig  8,  are  fitted  most  satis¬ 
factorily  by  A  =  0.8+0.05  rev/day  for  morning  conditions  (local  time  01-06  h) , 
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at  a  mean  height  of  200  km.  This  implies  an  east-to-west  wind  of  80  ±40  m/s. 

The  effect  of  meridional  winds,  p  rev/day,  is  expected  to  be  small,  as  the  orbital 
inclination  is  high;  however,  a  value  of  p  =  0.2  (south-to-north)  added  to  the 
A  =  0.8  curve,  has  been  calculated  and  plotted  in  Fig  8  (broken  line)  to  show 
its  effect.  The  average  meridional  wind  at  local  time  01-06  h  is  expected**^ 
to  be  about  40  m/s  south-to-north  in  Southern  latitudes,  equivalent  to  p  =  0.08; 
the  fit  is  not  improved  by  the  inclusion  of  p  .  Also  p  =  0.08  would  have 
about  the  same  effect  as  an  increase  of  0,03  in  A  .  Any  errors  due  to  neglect 
of  p  can  therefore  be  absorbed  in  the  error  of  0.05  in  A  . 

In  the  intermediate  months,  for  MJD  43282-43528,  the  values  of  i  could 
not  be  fitted  satisfactorily  by  a  mean  A-curve,  and  so  the  region  was  sub¬ 
divided  into  morning  (02-14  h)  and  evening  (14-02  h)  sections.  Although  the 
normal  time  division**  is  04-12  h  for  morning  and  16-24  h  for  evening,  it  was 
expedient  here  to  divide  the  day  into  two,  rather  than  four,  parts.  The  numeri¬ 
cal  integration  was  repeated  with  further  dates,  at  which  the  local  time  at 
perigee  was  02  h  and  14  h,  added  to  the  original  set  of  dates.  The  A-fittings 
were  made  to  each  section  separately,  including  inclination  'jumps'  at  14th- 
and  15th-order  resonances,  as  shown  in  Fig  7.  The  'morning'  sections  are  shown 
as  broken  lines,  the  'evening'  as  unbroken.  The  fittings  are  now  acceptable  and 
give  the  following  values  for  A  in  rev/day: 

(1)  MJD  43282-329,  A  =  0.9,  for  morning  conditions,  at  a  mean  height  of 

217  km,  for  May-June  1977 

(2)  MJD  43329-364,  A  =  1.2,  evening,  at  214  km,  for  July  1977 

(3)  MJD  43364-398,  A  =  0.9,  morning,  at  215  km,  for  August-September  1977 

(4)  MJD  43398-440,  A  =  1.2,  evening,  at  215  km,  for  September-October  1977 

(5)  MJD  43440-482,  A  =  0.7,  morning,  at  210  km,  for  October-November  1977 

(6)  MJD  43482-528,  A  =  0.8  ±0.1,  mean  (dotted),  at  217  km,  for 

December  1977-January  1978. 

The  sd  in  A  for  (1)  to  (5),  between  May  and  November  1977,  are  not  given,  as 

there  are  insufficient  inclination  values  to  determine  them  individually  (see 
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Fig  7).  However,  an  overall  sd  has  been  determined  using  the  THROE  program  , 
by  a  fitting  to  all  the  inclination  values  and  then  finding  the  change  in  A 
corresponding  to  a  change  in  the  di/dt  term.  The  sd  obtained  was  0.05. 

The  values  of  A  for  (1)  to  (5),  which  indicate  A  —  1.2  in  the  evening 
and  between  0.7  and  0.9  in  the  morning,  are  in  agreement  with  past  results 
obtained  from  other  satellites**;  however,  (6)  appears  to  give  an  anomalous  result 
for  A  .  To  date,  at  heights  below  400  km,  A  <  1  has  been  an  indication  of 
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morning  conditions.  For  a  height  near  220  km,  under  mean  conditions,  the 
expected  value  is  near  A  »  1.1,  te  a  wind  speed  of  about  40  m/s,  west-to-east. 

Here  (6)  gives  A  =  0.8  ±0.1  ,  which  represents  a  wind  speed  of  80±40m/s,  east- 
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to-west.  However,  this  apparently  anomalous  result  is  in  agreement  with  a  study 

where  for  similar  conditions  (height  215  km;  mean  latitude  57°S;  mid-summer; 

high  solar  activity),  the  wind  is  found  to  be  from  east-to-west  (Fig  6d  of  Ref  12), 

and  the  near  60  m/s  magnitude  is  in  close  agreement.  The  result  also  confirms 
.  1 3 

the  tentative  conclusion  from  analysis  of  the  Cosmos  462  orbit,  that  A  tends 
to  be  lower  in  summer  than  in  winter,  and  this  provides  a  basis  for  further 
refinement  in  terms  of  season,  in  future  analyses. 

7  CONCLUSIONS 

The  orbit  of  China  6  rocket,  1976-87B,  has  been  determined  at  51  epochs 
between  October  1976  and  February  1978,  including  orbits  determined  for  each  of 
the  last  16  days,  up  to  the  day  of  decay,  on  4  February.  The  total  number  of 
observations  available  was  more  than  4000,  and  there  were  Hewitt  camera 
observations  on  seven  of  the  51  orbits.  The  sd  in  inclination  varied  from 
0.0005  to  0.0028°,  while  the  sd  in  eccentricity  varied  between  6  x  10  ^ 

(equivalent  to  40  m  in  perigee  height)  and  3  x  10 

Values  of  density  scale  height  have  been  determined  from  the  change  in 
perigee  height  and  compared  with  values  from  the  CIRA  1972  reference  atmosphere; 
the  values  obtained  are  about  7  per  cent  higher  than  CIRA ,  on  average. 

Eight  values  for  the  atmospheric  rotation  rate  A  ,  in  rev/day,  have  been 
determined  from  the  variations  in  inclination,  for  heights  between  200  and 
230  km,  between  October  1976  and  February  1978.  Two  mean  values  (averaged  over 
morning  and  evening)  are  A  =  1.0  ±0.1  at  230  km  and  A  =  0.8  ±0.1  at  215  km. 

The  four  morning  values  are  A  =  0.8  ±0.05  at  200  km;  A  =  0.7  at  210  km;  A  =  0.9 
at  215  km,  twice.  The  two  evening  values  are  A  «  1.2  at  215  km,  twice.  The 
resonance  effects  on  inclination  at  14th-  and  15th~order  resonances  have  been 
accounted  for,  as  well  as  lunisolar  and  geopotential  perturbations. 

Seven  of  the  eight  values  of  A  agree  well  with  values  obtained  from 

other  satellites.  The  eighth  value,  A  =  0.8 ±0.1,  for  a  height  of  215  km, 

appears  to  be  anomalous:  it  is  equivalent  to  an  east-to-west  wind  of  80 ± 40  m/s 

whereas,  under  mean  conditions,  a  west-to-east  wind  of  40  m/s  was  expected.  The 

value  applies,  however,  for  December  1977-January  1978  at  a  mean  latitude  of 

57°S,  that  is  at  high  latitude  at  midsummer.  The  result  reinforces  a  previous 
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tentative  conclusion  that  A  is  lower  in  summer  than  m  winter  and  agrees  well 


With  the  numerical  model  of  Ref  12,  which  indicates  an  east-to-west  wind  of  about 
60  m/s,  in  similar  conditions  (ie  at  height  215  km;  mean  latitude  57°S;  mid¬ 
summer  and  high  solar  activity). 
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WL^bstract 

^"*The  orbit  of  China  6  rocket,  1976-87B,  has  been  determined  at  51  epochs  during 
its  17-month  life,  using  the  RAE  orbit  refinement  computer  program,  PROP  6,  with 
over  4000  radar  and  optical  observations  from  49  stations.  The  orbital  accuracy  is 
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The  rotation  rate  of  the  upper  atmosphere,  A  rev/ day,  for  the  height-band 
200-230  km,  was  calculated  from  the  decrease  in  orbital  inclination  (after  being 
cleared  of  perturbations)  to  give  th^ following  results: 

( 1)  cn for  morning  conditions,  ^  “  0 . 9  for  May-June  and  August-September  1977,  at 

215  km  mean  height;  (K  m  0.7  for  October-November  1977,  at  210  km; 

A  *  0.8  g£_p.05  for  January-February  1978,  at  200  km; 

(2) c^for  evening* conditions,  A  for  July  and  September-^ctojjer  1977,  at  215km; 

(3)  c^for  mean  (morning  plus  e'fefting)  conditions,  4*  1.0  . 1  between  October  1976 

and  May  1977,  at  230  km;  a  *0.8  $  0.1  for  December  1977-January  1978,  at 
215  km  and  mean  latitude  57*$. L  *  °r  ' 

Values  of  density  scale  height  have  been  obtained  from  the  variation  in  peri¬ 
gee  height,  including  several  values  during  the  final  16  days  to  decay.  Comparison 
with  CIRA  1972  values  shows  agreement  mostly  within  10  per  cent. 
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